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On the scattered field generated by a ball inhomogeneity of 
constant index in dimension three 

Yves Capdeboscq, George Leadbetter, and Andrew Parker 

Abstract. We consider the solution of a scalar Helmholtz equation where 
the potential (or index) takes two positive values, one inside a ball of radius 
e and another one outside. In this paper, we report that the results recently 
obtained in the two dimensional case in [Capl2] can be easily extended to 
three dimensions. In particular, we provide sharp estimates of the size of the 
scattered field caused by this ball inhomogeneity, for any frequencies and any 
contrast. We also provide a broadband estimate, that is, a uniform bound for 
the scattered field for any contrast, and any frequencies outside of a set which 
tends to zero with e. 



1. Introduction 

We consider a scalar field satisfying the Helmholtz equation with frequency 
lu > in R 3 . Given a prescribed incident field u l , a non-singular solution of 

(1.1) Au l + w 2 q a u z = in W 3 , 
we are interested in the solution u e € (R 3 ) of 

(1.2) Au £ + uj 2 q e u e = in M 3 , 

where, for |x| > e, u £ — u l + u s e , and q e equals q inside the inhomogeneity and go 
outside. We take the inhomogeneity to be a ball of radius e. The coordinate system 
is chosen so that the inhomogeneity is centered at the origin. In other words 

q if r < e 
go if r > e 

We assume that both g and g are real and positive. We assume that the scattered 
field satisfies the classical Silver-Miiller [MU169] outgoing radiation condition, given 

by 



1e{r) 



d (\ 

(1.3) — u% - iujy^u s e = o [ - 

where, as usual r :— \x 



or \r 
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The purpose of this paper is to provide sharp estimates for the scattered field 
u*, for any contrast q/qo and any frequency lj. The norms we use to describe the 
scattered field are the following. Given any / g C°(R 3 ), its restriction to the circle 
\x\ = R can be decomposed in terms of the spherical harmonics given by (|3.4[) . in 
the following way 



f(\x\=R) = J2 E fn,m{\x\ = R)Yl 
n— m— — n 

and f(\x\ = R) can be measured in terms of the following Sobolev norm 



(1.4) 



\f(\x\=R)\\ H - 



v E E l/n,m(N=ii)| 2 (2n+l)^, 
\ n=0 i 



) m— — n 



for any real parameter a. By density, this norm can be defined for less regular 
functions. For radius independent estimates, we shall use the norm 



(1.5) 



\ n=0 



It is easy to see that this norm is finite for a smooth / with bounded radial varia- 
tions. For a radial function, this is simply the sup norm for /. Finally, to document 
the sharpness of our estimates, we will provide lower bounds in terms of the semi- 
norms 



(!- 6 ) N p, 9 (/, K ) SU P SU P K 

p<n<q R>0 



: n+ 1 -(n+|) 5/ \2n+l) c 



, E \fn.m{\x\=R)\\ 
\ m— — n 



where q > p > are integers and k > 1 is a real parameter. These norms are satisfy 
the following inequality 

11/ (M = R)\\ H „ < Af a (f), and N- !00 (/, 1) < A^(J), 

and if for all R, f (\x\ = R) only has one non zero spherical harmonic coefficient, 



(/,l)=^(/) = sup ||/ (\x\=R)\\ Ha 



R>0 



They are the natural extension of the norms introduced in |Capl2] for the two 
dimensional companion problem. When the incident field is a plane wave, 

it z (x) = exp (i^/qouj^ ■ x) 

where £ is a unit vector in R 3 , for all R > 0, 



\u i (\x\=R)\\ H =^V) = 



^(2n + l)j n (uj E R/e f = 1, 

\ n=0 



whereas for any q > p and k > 1 , 

N; q (u\k) = max C{n)K n+ ^ n+ ^' 6 (2n + l) ff -* , 

where 2s > C(n) > 1 for all p and q (see (|3.30|) 1 . 

The motivation from this work comes from imaging. In electrostatics, the small 
volume asymptotic expansion for a diametrically bounded conductivity inclusion is 
now well established, and the first order expansion has been shown to be valid 



SIZE ESTIMATES 
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for any contrast |NV09j . It is natural to ask whether such expansions could also 
hold for non-zero frequencies, even in a simple case. Another inspiration for this 
work is recent results concerning the so-called cloaking-by-mapping method for the 
Helmholtz equation. In |KOVW10l |Ngul2^ ILZ11] , the authors show that cloaks 
can be constructed using lossy layers, and that non-lossy media could not be made 
invisible to some particular frequencies (the quasi- resonant frequencies). Within 
the range of non-lossy media, one can ask whether such 'cloak busting' frequencies 
are a significant phenomenon, that is, would appear with non-zero probability in 
any large frequency set, or on the contrary if they are contained in a set whose 
measure tends to zero with e. 

These questions were considered in two dimensions in |Capl2 . In this paper 



we show that these results extend, after some adjustments, to the three dimensional 
case. The proofs presented in this paper are very similar to the ones of the two- 
dimensional paper, but we believe the results, more than their derivation, could be 
of interest to researchers in various areas of mathematics. In numerical analysis they 
could be used as a validation test for broadband Helmholtz solvers, since we provide 
both upper and lower bounds for the scattering data. In the area of small volume 
expansion for arbitrary geometries, or in the mathematical developments related 
to cloaking, they provide a 'best case scenario' which can be used to document the 
sharpness of more general estimates. 

To make the results of this paper accessible to readers who are not familiar 
with Bessel functions, the main estimates are written in terms of the norms || • 
\\h" , N£ q (-, k) and A/" CT (-) introduced above, and powers of 2. Because no unknown 
constant C > appears in the results, this paper can be used as a black-box if 
the reader wishes to do so. Bessel functions do appear in one place, to describe 
quasi-resonances, but they turn out to be of the same nature as the two-dimensional 
ones, and so we refer to |Capl2 in that case. 



The main results of the paper are presented in Section [2j The proofs are given 
in Section [3] 

2. Main results 

To state our results, we introduce the rescaled non-dimensional frequency u> s , 
and the contrast factor A given by 

(2.1) u e := yfqooJE and A := */— • 

V 9o 

The following theorem provides our estimates for either small frequencies or for 
any frequency. 

Theorem 2.1. For any R > e, when max(A, l)u> £ < \ there holds 

(2.2) ||< (M = R)\\ H . < 2 4 / 3 |A - 1| w B I \\u* (\x\ = £ )|| ffCT _i . 

If Po > is the first integer such that the (po,m) spherical harmonic decomposition 
coefficient of u\ (\x\ — e) is non zero for some —po < m < po, then (|2.2[) holds for 
all max(A, l)ui £ < po + \ - Furthermore, for any R > e, when max(A, l)w e < i, the 
scattered field u s e also satisfies 

(2.3) \\u s s {\x\=R)\\ H „ 

< 2 2/3 |A - 1| max(A, !)"!;§( K(0)| + 2 2/3 max(A, l)w E W 7 ^^ (u l ) ) . 
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When R > max(A, l)e there holds 

(2.4) sup |K (\x\ = R)\\ H „ < 2 4 / 3 max(A, 1)75 AT (u 1 ) . 

Naturally, the variant of (|2.2p incorporating the more precise estimate given by 
(I2.3P for the Fourier coefficient corresponding to n — also holds by linearity. It is 
easy to verify that the dependence on uj £ in (|2.3[) is optimal by a Taylor expansion 
around e — (or uj — 0) for a incident wave u l with only one (or two) non-zero 
spherical harmonic coefficients for n — (and n = 1). Theorem 12.11 shows that this 
estimate is valid up to rescaled frequencies of order 1 when A < 1 , and of order 1 / A 
when A > 1. To prove the optimality of these ranges, we define, for t < 1, 

f 6 

n (t) = min neN such that t 2 < 1 - 



n>o\ ' ' (n+ ^) 2/3 

and for t > 1, 



ni(t) = min neN such that t < 1 



n>l( " (n+i) 2 /3 

Proposition 2.2. When A < 1 and e < R, 

sup ii< = i?)n ffCT > \ ^K:k™ («*. i) ■ 

When X > 1, R < eX for any integer q > ni(A), 

<"> „<L?<J K(H " H) »- s w"Ct (»' f »»G§ - 1 

In /act estimate (|2.5I) feo/ds i/ £/ie supremum is taken in the set 

0< Xu e <q+ 1.86g 1/3 + 1.04g" 1/3 . 

Note that ^ exp — l) > 1 when R < eX. 

To illustrate the sharp contrast between what these estimates show when A < 1 
and when 1 < A, let us consider the case of a plane wave. When A < 1, Theorem l2.1 
and Proposition 12.21 show that for any R > e there holds 

7 /9 * 4 < SU P K (W = R )Wh° < 2 4/3 ^ 
4 ^/2n (A) + l/i w>o -K 

When A > 1 , Theorem 12.11 shows that for any R > e , 

sup |K (|*| =£)|| ffl/3 <2 1 / 3£ 

and for any R> Xe, 



i>Aw e >0 " 



sup|K(N=i?)|| H1/3 <2 4 / 3 ^. 

The combination of these two estimates do not provide a bound when 2Aw e > 1 
for the near field e < R < Xe. Proposition 12 . 21 provides a lower bound in that case. 
For any p > n\(X) and any a, 

, (eX (R \\P+h-kP+h? l& 
jmpJK (M = B)\h- > 2-7/2 (^ ex P {Tx- 1 )) (2p + I)'"'. 



SIZE ESTIMATES 



5 



This lower bound grows geometrically with the upper bound of the interval of 
frequencies considered. In particular, for any A > 1, any R < Xe, and a 6 R, 

sup (\x\ = R)\\ H „ = +oo. 

This unbounded behavior of the scattered field is due to the existence of quasi- 
resonant frequencies, just as in the two-dimensional case. To characterize these 
quasi-resonances, Bessel functions are required. For t > 0, we denote by i? t (x) the 
Hankel function of the first kind of order t. The Bessel functions of the first and sec- 
ond kind of order t are given by J t (x) = 3? \ (x)j , and Y t (x) — 3 (x)^j . 

We denote by a u , m the m-th positive solution of J t (x) = 0. We denote by a^}} m 
the m-th positive solution of J' t {x) = 0. Finally, we write /3 Vj i the first positive 
solution of Y t (x) = 0. 

Definition 2.3. For any t > 0, the triplet (t, x, A) is called quasi-resonant if 

< x < /3 tA , 

and 

AJj (Xx) Y t (x) = Yl (Ax) J t (x) . 

The following proposition is proved in |Capl2| in the case when t is an integer, 
but the proof is unchanged for any t > \. 

Proposition 2.4. For any t>\ and X > at,\l fit.\, in every interval 

U t ,k = ^-^> ~A^ SUch that Ut > k C I 

there exists a unique frequency u>t,k such that the triplet (i,W(.fc, A) is quasi-resonant. 
There are no other quasi-resonances. In particular, no quasi-resonance exists in the 
interval (OjCt^/A), or when X < at^\/Pt,i- 

Since for any A there is only a countable number of quasi-resonant triplets, one 
could hope that outside security sets around the quasi-resonant frequencies, the 
scattered field could be bounded from above, even in the near field. This means 
excluding a countable union of intervals: a trade-off occurs between how much in 
the near-field one wishes to go, and how large the set of authorized frequencies is. 
The Theorem below is the result of such a trade-off. 

Theorem 2.5. For all X > 0, all e < 7~ 3/2 . R > e 1/3 , and a £ (0,1], there 
exists a set I depending on e, X and a such that 

\I\ <£ 1/3 |lne|, 

and 

Ifi f 1 / 3 

sup ||< [\x\ =R)\\ H „<- —N° +2+a {u 1 ) . 
v^we(o,oo)\/ " K 

If the contrast X is less than e~ 2 / 3 , this holds with 7 = 0. 
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3. Proofs of the main results 

Altogether, the conditions (|1.1|1.2I1.3[) imply that the incident field u\ the scat- 
tered field u| and the transmitted field u* = u e for r < e, admit series expansions 
in terms of special functions, namely 



oo n / \ 

(3.1) u\x) ~ E On,mJn(V®wN)Y™ f^J 

n — m— — n v 7 



(3.2) u»(a;) 



E 

n— m— — n 

oo 



) = T) VII/ 



(3.3) ^ E E On,mt»(w B ,A)j n (vM*IW ( u, 

n— m——n ^ '* 

In the above formulae, j n (x) — ^(hn^ {%)), and x — ¥ hn (x) is the spherical Hankel 
function of the first kind of order n, and Y™ are the spherical harmonics, given in 
terms of the polar coordinate 6 € [0, n] and (j> G [0, 2ir) by 



(3.4) Y™(M) = ^/^±l^_^|lp- (cos^exp(zm^), 

where P™ is the associated Legendre Polynomial. 

The reflection and transmission coefficients r n and t n are given by the trans- 
mission problem on the boundary of the inhomogeneity, that is, at r — e. They are 
the unique solutions of 

tn (w e , A) j n (Xu £ ) = j n (w e ) + r n (uj s , A) h$ (w e ) , 

Xt n (w e , A) j' n (Xw e ) = j' n (w e ) + r n (w E , A) h™' (w e ) , 

which are 

3? (h™ (w e ) Jn (\W S ) - Xj' n (XW S ) h£> (W s j) 

(3.5) r n {u e ,X) = ' ' 



hW'(w e )j n (XWe)-Xj n (XWe)hW(w e ) ' 

and, after a simplification using the Wronskian identity satisfied by j n (•) and 
(3.6) t n (w e ,X) ' ' 



w| feuv ( Us )j n (xw s ) - aj; (Aw e ) h<p (w £ ) ' 

In (13.11) . p.2[) and (|3.3I) . the ~ symbol is an equality if the right- hand-side is 
replaced by its real part, the fields being real. By a common abuse of notations, in 
what follows we will identify u % and u s £ with the full complex right-hand-side. 

To verify that this is the correct solution, we need to check that r n and t n are 
well defined. The fact that there is a unique solution to Problem (|1.2j) satisfying 
the radiation condition (11.31) is well known (see e.g. [Miil69 ). 

Lemma 3.1. h' n (w e )j n (Xw e ) — Xj' n (Xw e )h n (w e ) is non zero for all n € Z, u e > 
and A > 0. 

PROOF. Assume, for contradiction, j n (Xw s ) hn^' (uj e ) — Xj' n (Xw e ) hn^ (w e ) = 0. 
Then, as j n and j' n do not have common zeroes (see e.g. jOLBClOj ). either j n (Xw e ) 
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is non zero, in which case (|1.2p - (|1.3j) has the following solution 



(3.7) 



Ue (r):=i rofl^( A ^i) ifr<£ < 
U 11 ^) ifr>e, 



Or j' n (Xuj s ) is non zero, and (ll.2l) - (|1.3p has the following solution 
(3.8) 



uAr) . = J xtlfeM^f) if r < e > 
I W(u> c $ ifr>c. 



Both (|3.7I) and (|3.8[) would be solution of problem (11 .2|) - (| 1 . 3[) without an incident 
wave. There is of course another solution to that problem, U e = 0. Since Problem 
([L2 > ([L3 l) is well posed, see [Mul69 , we have a contradiction. □ 



We chose the three (semi-)norms M a and Np g (-,«;) because they are 

compatible with expansions (|3.1|) . (|3.2p and f|3 . 3[) . In particular, for any i? > 0, we 
have 
(3.9) 

1 

/ oc n n \ 2 

||«2(M=iJ)||^:=(X; E K, m | 2 (2«+1) 2CT r„(w e ,A)/iW ( We iJ/ e ) 
and 



vn-0 m— — n 



(3.10) Hu'di^eJll^: 
Whereas the other norms are 



E E Km| 2 (2n + l) 2<r |i„(u; £ 



vn— m——n 



(3.11) A/* CT := 

\n>0 m— — n 

and for any g > p > and n > 1, 



E E |a„, m | 2 sup|j„ (x)| 2 (2n+ F 

z — ' — ' x>0 



2a 



(3.12) N£, (u\/e) := sup 

p<n<g 



+|-(«+|) 5/6 ( 2n+ l)- SU p| Jn ( x ) 



a;>0 



\ 



m——n 



3.1. Proof of Theorem 12.11 The two dimensional results found in |Capl2| 
are easily translated into three dimensional ones for the following reason. 

Proposition 3.2. Let R n [ui e , A) be the reflection coefficient associated to prob- 
lem (|1.2[) posed in a disk of radius e in dimension 2, with the appropriate out-going 
radiation condition. This reflection coefficient is defined (by the same formula) 
when n is an arbitrary positive real number, and uj £ is replaced by any real x > 0. 
Then, for any n > and any x > there holds 

r n (x,X) = R n+ i(x, A). 

Notation. From this point onwards, we use the short-hand v to represent the 
number n + i. 



Proof. For n > we introduce 



j v {x) = - , and k v {x) = 



v J v (x) ' 



v Y v (x) ' 
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Then, the reflection function R v introduced in two dimensional problem considered 
in |Capl2| is 

R (x X) — g v (Xx)-g v (x) 

g u (Xx) — g v (x) + itaxi6 1/ (x)(g v (Xx) + k v (x))' 

where i 2 = — 1. For any v > and any x such that J u (x) ^ 0, we write 

tan(Ux)) = ^M. 

J v (x) 

The properties of the function R t were studied in |Capl2| for any t > 1. When t is 
an integer, R t is the reflection coefficient associated to problem (ll.2j) in dimension 
2. Note that the identities 

(3.13) j t _i (x) = y^Jt (x) and y t _i (x) = \J^Y t (x) , 

valid for any t £ R yield that 

tan(^(*)) = ^4. 

Jn (x) 

We introduce the notations 

u n = Xj' n (x) jn (Xx) - \xj' n (Xx) j n (x) 

v n = xy' n (x) j n (Xx) - Xxj' n (Xx) y n (x) , 
Then, the reflection coefficient corresponding to (|1.2I) is given by 

r n (x, X) = . 

If we introduce for any t € K, outside the zeroes of j t _i (x) and Vt-± (x), 

x J't-i W x y't-i- ( x ) 

ft(x) = - 2 . r , and s t (a!) = 2 —-^, 

tj t -i(x) ty t _i(x) 

we can rewrite, when Xx is not a zero of j n jn (Xx) or j/ n (x), 
u n = -nj„ (a:) j„ (Ax) (/„(Ax) - /„(x)), 
v n = -ntan0„(x)j„ (x) j„ (Ax) (f v (Xx) + s v (x)). 

From (13.13)) it follows that for any t > 0, 

ft(x) = gt(x) - — and s t (x) = k t (x) + 
where these functions are defined. This in turns implies that 

U n fv{Xx) - f v (x) 



f v (Xx) - f u (x) + itan0„(x) (U(Xx) + a v (x)) 
g v (Xx) - g v (x) 



g v (\x) — g v (x) + itan0„(x) (g v (\x) + k u (x)) 
= R v (x,X), 

and this identity extends to the zeroes of j n (x), j n (Ax) or y n (x) by continuity. □ 
The following Lemma then follows mostly from |Capl2| . 
Lemma 3.3. For any A > and n > 0, 
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(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



For x € (0, min(a[ / 1 ]/A, fiv,i)) there holds 



2 3 Jn (#) • 



For a: G (0, min(l/ A, we have 



r n (x,X)hW(x) < 2*|A - lj 



(2l/)V3 



jn (x) . 



Furthermore, when < max(A, l)x < h 



rQ{x,X)h^ ] (x) < 23|A- 1| max(A, l)s 2 , 



and w/ien < max(A, l)a; < ^ and n > 1, 



r n (a:,A)/iW < 2 I |A - 1| max(A, if 



When A < 1 and n > rto(A), we have 



(2i/) 1 /3 



Jn (1) ■ 



• When A > ni(A) iftene eiisis x n < a^i/A suc/i tta< 

(3.19) |r„(a;„,A)| =1. 

Proof. The existence of x n satisfying (|3.19[) follows from Proposition 
Thanks to Proposition 13.21 and because for any n > 0, we have 

jn (x) J v (x) 



h^ (x) Hi 1} (x) ' 

the inequalities (13. 14|) , (|3. 15|) and (|3.18[) are proved when n > 1 in |Capl2| , Lemma 
7.1. We will now check that (|3.14j) holds when n — 0. From (|3.5j) there holds 



Since 



ki,A)|<l. 



Ji (x) = — t&n(x)Yi (x) , 



we have for all x € (3/5,7r/2) 



(3.20) 



r (x,X)h^ (x) <] n (x) ( 1 + 



tan(3/5) 



The zero is given the first positive solution of tan(x) = 2x. It satisfies 1 < 

2 i 1 

oA 1 ^ < /3i ;1 = |. We can thus consider only < x < min(|, a^A -1 ). 
We have from Proposition [3T2l for n = 0, 1, 2 

r„(x,A)/iW (x) 

(ff n+ 1 (Ax) - 5„ + 1 (»)) (1 + i tan 6» n+ 1 (a:)) 



= ~jn (x) 



3„+i(Ax) -g n+ i(x) + itan0 n+ i(ic) (g n+ t(Xx) + k n+ i(x) 

jn (X) S„(x), 
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and using the Wronskian identity satisfied by J n+ i (x), Y n+ i (x) and the recurrence 
relations satisfied by Bessel functions, we obtain (see |Capl2| for details) that 

_ U n (x) 
s n\ x ) / \ . i 

U n (X) + 1 

with 

u n (x) = [2n + 1)jH^i (x) J n+ i (x) (g n+ i(x) ~ c? n+ i(Ax)) . 

Note that g n+ i(x), J n +± ( x ) an d l-^+i ( x ) I have simple analytic formulae. For 
example 

X . II. 



In particular, it is easy to show that <7„ + i is decreasing on (0, o&_i x ) 



When n = 0, we have the following bound for all < x < min(3/5, ai A 1 ). 

2 ' 

|«o(af)|< max |jlfi 1) (x)Ji(x)| max \gi(x)\ < -. 



4,i 



As a consequence, 



i i \\ s' KO) s , 

Sola; < -Tj i 7 777 < 1- 

i - K x) 



Which concludes the proof of Q3.14p when n = 0. 

For any n = 0, 1, 2, and max(A, l)x < n + |, we have 

|u„(x)| < max | 2n 1 W H„li (x) J n+ i (x) [ max \g n+ i(x)\ 

2n + 1 , r 
= — 2~ M n <l, 

where M„ = max 0<J . <n+ i \(2n + l)f (x) J„ + i (x) | for n = 0, 1, 2. Therefore 

(3.21) | S „(x)| < < 2 Mx)|. 

|l-|u„(x)|| 2-(2n+l)M„ 

,2 i 

To proceed, note that -j-^gi is negative and decreasing, and for x < n + 



2 



d 2 fn + l 



Thus a Taylor expansion shows that for < x < n + | , 

(3.22) |w«(x)| < M„k„|A - 1| max(A, l)x 2 . 

Inserting (I3.22[) in (13.21)) together with the values of Mq and kq, we obtain 

|soO)| < 23 |A- l|max(A, l)x 2 

for all < 2max(A, l)x < 1, which implies (|3.15l) . and (13.161) since \jo (x) | < 1. 

Inserting p. 221) in p. 211) together with the values of M n and K n for i = n and 
n = 2, we obtain for x < n + ~ 

(3.23) |s„(x)| < \/3|A- l|max(A,l)x 2 (n+ -J 



SIZE ESTIMATES 



11 



For all n > 0, it is known that Par84 for < x < y < a„\ 



(3.24) 



jn (x) < —in (y) exp 
y 



2 9 

y — x 
2n + 5 



In particular, for all x < 1, we have 



(3.25) ji (x) < xexp - j 1 (1) , and j n (x) <x exp - j n (1) for n > 2 



7 

Combining p~23|) . (1X2"5]) and (j3~T5|) we obtain ((3~T7|) . 
We may now conclude the proof of Theorem 12. II 



□ 



Proof of Theorem 12.11 For convenience we write A := (\x\ = R)\\ H a 
Formula (13.91) shows that 



oo n 



[0J E 



n—po m——n 



h { n ] (u e R/e) 



Note that a; ftn (a;) is decreasing (see e.g. }Wat95| § 13.74), therefore for all 

R> E, 

2 oo n 2 

(3-26) A 2 <^Y. E |«n,™(2^) c7 r„H,A)^ 1 )H) . 



n— po in — — n 



Estimate ()3.15jl in Lemma [3T31 shows that when max(A, 1)oj £ < po + i < ia 

(2i/)" |r„ (w 8 , A) (Wb )| < 2 4 /3| A „ i| We ( 2l/ )'-i |j„ ( We )| . 

Inserting this bound in (|3.26l) we obtain 

A<2 4 / 3 |A-lK^||^(N=e)||^_i , 

which is estimate (|2.2[) . Using p,17[) instead, we obtain for n > 1 and u; e < 1, 

(2z,) 2 - |r„ ( Wej A) (o; e )| < 2 4 / 3 |A - 1| max(A, 1) 2 ^V^ |j„ (1)| 2 . 

Inserting this bound in p. 261) together with (13. 16)) we obtain 

A < 2 2 / 3 u 2 \X- l|max(A,l)-| (|a | + 2 2 / 3 w e max(A, lJAf* - * (u 1 )) , 

which proves (|2.3[) since |ao| = |u J (0)|. Let us now turn to (12. 4[) . Formula p. 9 
shows that 

oo n 2 , (1) / t> / \ 

A2 =H E «n, m (2^rr n K-,A)/ l ( 1 )(max(A,lK) "" ' " " " 1 



n—po m——n 



hn^ (max(A, l)to e ) 



When R > max(A, l)e 7 since x \ hn^ (x)\ is decreasing 



oo n 



(3.27) A 2 < (max(A,l)-|) ^ ^ |a„, m (2z/) CT r n (w e , A) (max(A, l)w e ) 



n—po m=—n 



We define two sets of indices, 

I := |n e N such that w £ < min(a^iA -1 , /3 u ,i)\ 
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and J = N \ I. S incG fin (x) is decreasing, estimate (|3.14l) shows that for all 
n £ I, we have 

r„ (w e , A) /#> (max(A, l)w e ) < r n (u e , A) (u> e ) < 2 4 / 3 sup \j n (x)\ . 

x>0 

On the other hand, v < a„ \ < /3 Vj i, therefore when n 6 J, max(A, l)(J e > o^\, 

' fll 2 

and it is known (see [Wat 95] § 13.74) that when x > n. ^x 2 -n 2 Im 1 ' (x)\ is an 

increasing function of x with limit 2/tt. Furthermore, it is also known (sec [Wat95 
§ 15.3) that for all n > 

Therefore, since \r n (u) e , A) | < 1, we have 

r n (u E , A) hW (max(A, l)w e ) ' < (max(A, l)u e ) 



On the other hand 
2 



7rmax(A, l)w e 



ff« (max(A, l)w e ) 



< 



7r max(A, l)w £ 
2 



ffW (max(A, l) We ) 



J 1 ) 



< 



2 3 1 



^;i7rV(«S) 2 -^ 

As it is known (see [LanOOj ) that for all n > 
(3.28) 



2 7TQ! 



(1) 5 V 2 / 3 ' 



I'll 



0.539 < v 1/3 sup J„ (x) = v 1/3 J„ (ctfl ) < 0.675, 

x>0 ^ ' ' 



this shows that 
2 



7r max(A, l)w e 
We have obtained that, for n G I U J, 



H« (max(A,lK)| 2 < — ^J,, (ag) 2 - |i„ (og) 



(3.29) 



»"« (w e , A) /#> (max(A, l)w e ) < 2 4 / 3 sup |j„ (x)| 



z>0 



□ 



Inserting this bound in (|3.27|) . we obtain (|2.4|l . 

3.2. Proof of Proposition I2T21 To prove Proposition 12. 2[ we shall use the 
following intermediate result, related to Bessel functions. 



Lemma 3.4. For any n > 0, there holds 
(3.30) 



n ■ i sup \j n (x) I > 5/6 > 0.663 sup \j n (x) | 
2sm~x>o (2v) b ' h x>0 



For any n > no (A) there holds 



(3.31) 



R 



> 



e 0.58 



fli/Ve j„M x>o 



sup|j„ (at) |. 



For any A > 1, n > 1, and < 1, £/iere Zioids 



(3.32) 



) ( ^ a ".i 



> 



2 -7/2^_/fA 

(2^)i/3 exp veA 



.5/6 



SUp|i„ (x) |. 

x>a 
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We prove this lemma below. We can now conclude the proof of Proposition ^. 21 

Proof of Proposition 12.21 Let us start with the case A < 1. Starting as 
before from formula (I3.9[) . we have 

sup \\u 8 s {\x\ = R)\\ Ha 



ui e >0 



> sup sup V] \/|<V m | 2 r n (w e ,A)/i£ 1) (w e ) (2z/) cr 

w s >0„>„ (A) m= _ n V 

{vR/e) 



h { n ] (ujsR/s) 



h£ ] (w e ) 



> sup y|ffln,ml 
n>np (A) 



ytrn|^n(^ A)/li 1) (l/)| 2 (2l/) e 



Using now the bounds (|3.31|l in Lemma 13.41 and (|3. 18|> in Lemma 13. 3[ we have 



{2v)° r n (v t \)hW(v) 



{vR/e) 



1 £ l 

> 7 p(2fr~*sup lin(aj)|, 



Therefore 



1 e 



sup |K = > ~ sup /|a„, m | 2 sup |j„ (fl)| 2 (2i/)*-*, 

w E >0 4 - tt -n>n (A)V - R>0 

as claimed. We now turn to the case A > 1. Note that for all q > 1, using the 
upper bound given in (|3.34l) we see that a q+ i x < 6q. Therefore, choosing for each 
q the frequency x q given by Lemma 13.31 we have 



sup =R)\\ H „ > sup 

0<Aw e <6g K 




R 

Xq 

£ 



■ , \a n ,r. 

-1 \ m=-n 



{2v)° 



A e 



{2vy. 



□ 



The conclusion follows from estimate (|3.32j) . 

Proof of Lemma 13.41 To estimate the maximum of \j n (x) |, we proceed as 
follows. Note that the maximal value occurs at the first positive solution of j' n (x) — 
0, which we will denote j v> i. We compute that 

• / / \ Jy (x) f xj' v (x) _ 1 

Jn[X) ~ V 2 a;3/2 ^ J v {x) 2 / 
It is known (see e.g. |Capl2| Proposition A.l) that for x < v and v > 1, 

JJM > *!L > 1 

J„(x) V2 2' 
Therefore 7„ i > v. This implies that 

2n x> o \ In V / x >o 

Next, note that from (13.28p . and the bound v < OLy x we have 



> 



- i/3 j, («s) 



0.675(2,) 5 /a - .6752 i/a^V 2n > °- 663 ^ ljn {x) 



> 
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On the other hand, v — > v 1 ' 3 ' J u (v) is an increasing function (see |Wat95j § 15.8), 
therefore 



2 sin — = \priJi 
2 i 



(I) < sfrJv M (21/) 1/3 = (2z/) 5 / 6 |j n (i/) | < (2i/) B /« sup | in ( X ) |. 

It is well known (see e.g. |OLBC10] ) that for all x > 0, (x) | > A. It is also 
known (see |Wat95j § 15.8) that v — > — j^j^y is a decreasing function. Note if 
n > n (A) then z/ > 6 3/2 . Therefore 

M I < in M J i + ^77^ < 2 - 01 ^ M • 

V Jio(lO) 

(u) 



Combining these two bounds we obtain that 



(3.33) 



R 



> 



e 1 



#2.01*/ j n (u) 



Together with f|3 . 33[) this shows that 



hi 1 



R 



> 



e 0.58 



D 1/6 • ^ SUp|j„(x)|. 

Let us now turn to (|3. 321) . It is known (see |QW99| ) that 

3 

ib 

Where A\ is a universal constant, s» 1.855757082. 



(3.34) 



oc v ,i =v + c(^ 1/3 where A Y < c(v) < Ai + —A\v 2/3 . 



Let us first assume that j| ( 1 + M ( ^ 



c(u) \ c(v) 



< 1, where M is given by 



(3.35) M(x) = 2 + 2x + x 5/4 . 

Then, using p.34j) we find 

R l + c{v)v- 2 ' 3 

It is known (see )Par84j ) that for all x < v, 

J, (x) < g exp (l _ J„ („) < _^ (5 exp ( 



1 



Therefore 



where 



In i 



■I, ( ! ■ 



R 



1 



2i/V3 1 eA 



1? 



1? 



— exp 1 - 



eX 



,,6/6 



eA 



I? 



eA 



= ^ ( £ ) + - 5/6 ( 1 - £ ) + - ^ ( ^ ) + - «M) 



i? 



eA 
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Since x — > (1 — x) 1 ln(a;) is increasing for x < 1, we have 



R 



eA 



1 - - In ^ 



1? 



eA 



c(i>) 



< 



In 1 



c(y) 



Using the definition of M (|3.35l) , we see that the right hand side of this last inequal- 
ity is an explicit function of v~ 1 > 2 c(j/), which is negative when z/ _1 ' 2 c(i/) < 2.18. 
Using (ETM|1 , we see that for v~ x l 2 c{y) < 2.16 for all v > 3/2. Thus for all v > 3/2, 
there holds 



■/,- ( ^a„,i j ; 



1 



2i/V3 VeA 



i? 



i? 



— exp 1 - 



eA 



Next, we note that x — > —J v (x) Y v (x) is minimal on (0, v) at x — for all v > 1, 
where it equals (7tj/) _1 . Therefore 



eA ^j* , 



,,5/6 



and, using that 



I —a i 



R 



i? 



we obtain 
(3.36) 



ft! 1 



eA 



> 



2 1 



7T ^ 7 /6 I i? 



eA 



■ exp 



R_ 

7x 



eA 



5/6 



Let us now suppose that 1 > > ( 1 + M 



Xe 



Since x 



x 1 exp (2 — 1) is decreasing when x < 1, we obtain an upper bound on 



eA 



R 



R 



— exp — - — 1 



■A 



by replacing R/(eX) by its lower bound and civ) by its upper bound, given in 
()3.34[) . The resulting expression is an explicit increasing function of n, with limit 
exp(2A\) « 980. It is then possible to verify by inspection on a finite range for v 



that for all v > |, 



(3.37) 



/V3 \, i? 



eA 



R 



— exp — - — 1 



eA 



,V3 



< 11, 



this inequality being automatically satisfied when v > 10 6 for example. The max- 

JL 

sX 



imum occurs near v — 2838. On the other hand, ^(%v,i < cvi and x — > /i^ (a;) 
is decreasing, and see e.g. |OLBC10] 



-y v K,0 > ^^ 1/3 , 
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therefore when x < a Vf i 
(3.38) 



hi 1] (x) = Jtt—Y v (a Vl i) > L , 1 = 



> 



2n 5/6 ■ 



The bounds (|3"757|l and (13381) show that when 1 > -§ > (l + M (jt^ , 
we have 

R \ 1 1 /eA 



(3.39) 



cA 



> 



15 z/ 7 /6 I i? 



- exp 



,,5/es 



Combining p.36j) and (|3.39[) we have obtained that (|3.39|) holds for all A > 1, all 
n > 1 and all R < eX. To conclude, note that using (|3.30|) we have 

77 "T76 > 2 ' 7/ %o u/3 sup l J " ^ ■ 
15 ^ 7 / b (2^) 1 /J x>0 

□ 

3.3. Proof of Theorem [2751 The proof of Theorem [23] follows the line of the 
proof of the corresponding result in the two-dimensional case proved in |Capl2| . 
The first step is the following proposition. 

Proposition 3.5. For any < r < j and A > 7, we define 

< ^- SUp \j n (x) 
6 T x>0 



then 



A n (T~) — |cj e > such that r n (cj S) A) (u £ ) 

B n {r) = (0,oo)\ A n (r), 
B n {r) C |J I Vtk (r), 

keK(X,n) 



where /^(t) is defined by 

(3.40) Iv,k( T ) '■= i x e Uv,k suc h that \g u (Xx) + k v (x)\ < r |fc„(a;)|} , 

and where K(X, n) is the set of all positive n such that a", ^ < v\. Furthermore, 



(3.41) 



, 2vhxX 
\B n (r)\ < 4t^—. 



When n — 0, the same result holds for r < | . 

PROOF. First, note that Lemma [3731 shows that (0, a„\/X) c ^,fc(T). 

Furthermore, we have shown in the proof of Theorem 12.11 that when uj £ > a^\, 

r n (u e , A) hn^ (u) e ) < 5/2sup :E>0 \j n (x) |, thus (a^,oo) C A v ,k(r). 

Thirdly, using the bound \r n (u) e , A) | < 1, we see that when v < uj s < ct^\, we 
have 



r n (u) e , X) hn' {w e ) hfi J (ui e ) 
< 



su P a; >0 \jn (X) | 2 \j n (We) 



Y v (uj e ) 2 
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( 1 ) 2 2 2 2 

On z/ < u) e < a u {, Y v (ui e ) jJ v {io e ) is decreasing, and Y v (v) /J v {v) is a decreas- 
ing function of v. Therefore 



2 



, <1 + -4ttV<5, 



SUP ;C>0 \jn (X) | 2 4(f)' 

and we have obtained that 

Bnir) C (<£}/\,v) . 

Next, it is known (see |Capl2| , Proposition 8.3) that when uj e € (0, v) and u £ £ 
Iv,k{T~) for some fc, then x € A n (r) when n > 1. 

The argument is simple. It turn out that by a simple calculus argument using 
the formula for r n , when g v (x) > 0, fcj,(x) > 0, and k v {x) > ^g v {x), \g u (Xx) + 
k v {x)\ > Tk„(x) 7 then x € A„(t). When a; € (0,^), g„(x) > 0, fc„(x) > 0, and 
k w (x) > ^g u (x) therefore the inclusion holds. Since gi(x) > 0, ki(x) > 0, and 
ki (x) > gi (x) on (0, |), the same is true when n = 0. 

The proof of Proposition ^. 5l will be complete once estimate p. 411) is established, 
for A > 7. Since it is proved in |Capl2 Proposition 8.2 for n > 1, we only need to 
consider the case n = 0, and w e G (o^/A, |). We have 

2x 

01 fx) = H : — and ki (x) = 1 + 2a; tan x > 1 + 2x 2 , 

2 V 7 tanx ^ 1 

and a.i fc = kix. Introducing 

(3.42) O (0,f)\U fe {WA} -> R 

31 (Ax) 

x -> 



fci(x) ' 

2 v ' 

we have 

0o(/i lfc (T)) = [-l-r,-l + r]. 

We first verify that <pi is one-to-one on Ii fe (r), for r small enough and A large 
enough. Differentiating <fr we find 

2t „ oi(Ax)+fci(x) / 1 — 01 (Ax)fci (x) \ 

= ir§> (1 - A > + ' j ( x)' ( ' „ ' - 2 *J • 

When —gi (Ax) > fci (x), we have 

1 -, 9 i(Ax)fci(x) 2 + 4x 4 + 4x 2 n l + 2x 4 n 
2— 2x > 2x = — > 0, 

therefore 

2 v ' 

When — gi (Ax) < fci (x) and x € Ii fc(x), we have < (71 (Ax) + fci (x) < rfci (x), 
and 

.91 (Ax) + fci (x) /l- 5 i(Ax)fci(x) \ T l-4x 2 + fc|(x) &Tx) 2 

2x < — — — < 



fcf(x) I 2x J xfci(x) 2 5 fci(x)' 
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since 1 < < Xx. Finally, A 2 < 1.05(A 2 - 1) when A > 7, thus for any r < 3/4, 

2 I 1 

we have obtained that 

^)<-^(l-A 2 )<0, 

2 v ' 

for all x £ Ii j.(t). In particular, if Ii k {r) — [Cfcj^fc], we have 

2r= ~cb'(t)dt>(\ 2 -l) [" -L^dt = (X 2 -l)\I hk (r)\-^— P -r^dt. 
Therefore 

/ c ^ - 1 



di 



fee.ff(A,o) feeA'(A,o) V* 5 V y 

since x/k{x) is increasing on (0, 1/2), and ?7fc > cSP,/X, 

2 ' 

E ^-^fr^) * E fc i(^ri a) 

fceA(A,0) V" 7 * 2 w / k£K(\,0) \ / a i,fc 

a n i 

< max-TTT ttt / —ki(t)dt, 

where oS^\ = - — Wtt is a convenient choice (but any number smaller than a , , 
and greater than zero would do to write the Riemann sum) . The distance between 
two distinct positive solutions of tan(x) = 2x is at least it, therefore 

A A A 2 -l 

maX 7TT 7TT = TTT 7TT < 2 



fcfA (!) (!) ~" (!) (1) A ' 



since A > 7, and 



1 a (i) 

2 -k h {t)dt = ln(A) + 21n(cos(-|^)) < In(A). 

i» 

We have obtained that 

E. _ , ,, At In A 
|/ ilfc (r)| < -jj-, 

fcGA(A,0) 

and this concludes our proof. □ 



The second step is to use Proposition ^. 5l to derive an estimate for the scattered 
field, for large enough contrast. 

Lemma 3.6. Suppose A > 7. Let ry max be the following decreasing function of 
the contrast 

5 In A 

(3-43) r?max = - — . 

Given a > 0, for any rj > such that 

1 

V < -?7max 

a 
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there exists a set I depending on -q, a, e and A such that 

V 



\I\ < 



and, for any R > e 



sup |K (M = R)\\ H . < 18 | *EZ.Ar°+*+° („<) . 
Proof. Recall that we established in (|3.26p . that for all R > £ we have 



oo n 



k cm = #)ii^ < ^ E E ( 2 ^ r « A ) ^ ) ( w =) 

n—po m— — n 

For < a < 1, and r/ < ~T7 roax , let t„ be given by 

1 m/ 
T " = 4(2^) 2 +^ max - 
If cj £ ^ U^Li Bn,r n , we have thanks to Proposition 

oo n 



E \an, m (2v) a r n (uj E ,\)hM(cj e ) 



n=po m——n 
oo n 

< E E 



n=0 m=—n 
//fl- 



an,™ ^Z^)" 7 - Sup \j n (x) 

2 T n X >Q 



A^ +2+Q K) 2 . 



From Proposition 13. 5[ we also know that 



OO 9 

a -k 

A » A r? ma x^;(2n + l) T ^ " 77 10 " ^ 

n— 1 n— n— v 7 



U° „ 4 In A v-^ , * In A ?7 v-— 

s Bj t„i < E( 2iy ) T » = i — - E 



To conclude, note that the set of excluded frequencies for ^/qou> is | U^Li ^n,r„ ■ D 
We can now conclude the proof of Theorem [ 



Proof of Theorem 12.51 When < A < 1, Theorem 00 implies Theo- 
rem EL3 with 7 = 0. When 1 < A < £~ 2 / 3 . Theorem O shows that for all 
R > e 1 / 3 , we have 

su P K(|.*|=i?)||^ <2^ £ -^-N° 

so we can again select 7 = 0. Suppose now A = e s , with s > 2/3. Then A > 7, and 
we can apply Lemma 13.61 Choosing 

3 



we have for all < a < 1, 



77 = - se s+3 |lne| 



1 

7? < -l)mnx, 



and there exists a set 7 depending on A, a and e such that for any R> e, 
su.» ||< (|x| - 7?)|| ff „ < — A^+ 2+Q 

v /^u£(0,oo)\/ a " 
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The size of the set I is bounded by 

< |se s - 1/3 |lnej < e 1 / 3 |lne| , 
since s — > se s_1 / 3 is decreasing when s > 2/3 > |ln£r|~ . 
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